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Abstract: Stochastic quantization can potentially be used to simulate theories with a 
complex action due to a nonzero chemical potential. We study complex Langevin dynamics 
in the relativistic Bose gas analytically, using a mean field approximation. We concentrate 
on the region with a Silver Blaze problem and discuss convergence, stability, fixed points, 
and the severeness of the sign problem. The real distribution satisfying the extended 
Fokker-Planck equation is constructed and its nonlocal form is explained. Finally, we 
compare the mean field results in finite volume with the numerical data presented in Ref. 
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1. Introduction 



Theories with a complex action are not easy to solve numerically since approaches based 
on importance sampling break down. This is commonly referred to as the sign problem. An 
important theory in this class is QCD at finite baryon chemical potential, with a complex 
fermion determinant satisfying [det M(|u)]* = det M{—fi).^ Several methods have been 
devised to circumvent the sign problem in QCD, mostly at small chemical potential and 
in the vicinity of the crossover between the confined and the deconfined phase |2|, ^, ^, |5|, 
H) 0) @) @) ^ 0' HI' HI' US' -^^^ ^ detailed lattice QCD study of the sign problem at 
small chemical potential, see Ref. |]l5|. Considerable insight in the QCD sign problem has 



also been obtained with Random Matrix Theory [|T6|, ||, |T|, |T|, |o[ H . In some theories 
the sign problem can be eliminated altogether, using a reformulation in terms of different 
degrees of freedom 

does not rely on importance sampling, it can po- 



Since stochastic quantization 
tentially be applied to theories with a complex action using complex Langevin dynamics 



Studies in the 80's, however, have given mixed results, see e.g. Refs. ||28| , 



For an extensive review and more references, see Ref. |30|. Recently the approach was 
reconsidered as a method to solve nonequilibrium quantum fields dynamics in Minkowski 
spacetime |31, 32, 33|. It was shown that instabilities, which plagued earlier studies, can 



^In case of a complex chemical potential, this relation becomes [detM(/i)]* = det A'f(— /i* 
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be controlled by using small enough Langevin stepsizes. Moreover, insight in the conver- 
gence properties of the method can be obtained from features of classical flow diagrams. 



Other recent applications include PT symmetric theories |34] and unbounded actions p5| . 
In Ref. |^6[ we applied stochastic quantization to various theories with a nonzero chemi- 
cal potential. In particular, we considered QCD with static quarks, in which the fermion 
determinant is approximated but the full gauge dynamics is preserved. First results on a 
4^ lattice are encouraging. The required extension from SU(3) to SL(3,C) is discussed in 
detail. 

The sign problem in QCD at finite chemical potential does not arise because of the 
anticommuting nature of the quark fields. Also in bosonic theories with a nonzero chemical 
potential and an action that behaves under complex conjugation as S*{^) = S{—fj.), the 
sign problem appears. In Ref. |l| we considered the relativistic Bose gas (a self interacting 
complex scalar field) in four dimensions in the presence of a chemical potential as one of 
the simplest examples of a relativistic field theory with a severe sign problem. Like QCD, 
this theory has a Silver Blaze problem : at strictly zero temperature and small chemical 
potential, bulk physical observables are independent of the chemical potential, even though 
it enters explicitly in the microscopic dynamics. At larger chemical potential, the system 
enters a Bose condensed phase. The fj, independence below onset and the formation of a 
state with nonzero density above onset is similar to what is expected to occur in QCD 
at zero temperature. It was demonstrated in Ref. Q that complex Langevin dynamics 
reproduces the expected physics, on lattices of size TV^, with N = 4,6,8,10. The sign 
problem was shown to be severe. However, no obstacles related to the sign problem, the 
Silver Blaze problem, or in taking the thermodynamic limit were encountered. 

In this paper we complement the numerical study of Ref. with a detailed analytical 
study in the mean field approximation. We concentrate on the region with the Silver Blaze 
problem. The paper is organized as follows. In Sec. ^ we remind the reader of the model 
and the corresponding complex Langevin equations. In order to prepare for the mean field 
analysis, we first discuss the case without interactions. In Sec. ^ we summarize the exact 
results in the free field limit, using standard field theory. Subsequently the free Langevin 
equations are solved analytically, both for continuous and discretized dynamics. We discuss 
convergence and stability properties. In Sec. § the stationary solution of the Fokker-Planck 
equation is given, again ignoring interactions, and shown to be in agreement with the 
solution of the Langevin equations in the limit of large Langevin time.^ The criterium 
for convergence is derived from this distribution as well. We find that the real probability 
distribution is highly nonlocal and explain why. In Sec. ^ interactions are included and the 
analysis is extended to the mean field approximation. We derive fixed points of the mean 
field Langevin equations at finite Langevin stepsize. Finally the mean field predictions in 
finite volume are compared with the nonperturbative results obtained by complex Langevin 
simulations M. The appendix contains a short remark about lattice dispersion relations. 



^Recently, an interesting approach to study stationary solutions of complex Langevin dynamics was 
presented in Ref. |3g]. 
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2. Relativistic Bose gas and Langevin dynamics 



We consider a self-interacting complex scalar field in the presence of a chemical potential 
with the continuum action 



5 



+ /i(</.*a4(/.-a4</)*</.) + A|0|'^ 



(2.1) 



The euclidean action is complex and satisfies 5'*(/u) = S{—fi). We take > 0, so that at 
vanishing and small the theory is in its symmetric phase. 
We study this theory on the lattice, with the action 



{2d + m^) (j)l4ix + A I 



X+V + 



(2.2) 



As always, chemical potential is introduced as an imaginary constant vector potential in 
the temporal direction |^9|. The number of euclidean dimensions is = 4, the lattice 
spacing aiat = 1, and the lattice four-volume is 17 = N^Nr, where Ng (N-r) are the number 
of sites in a spatial (temporal) direction. We use periodic boundary conditions. 

In order to formulate the complex Langevin equations for this theory, the complex field 
is first written in terms of two real fields (pa (« = 1,2) as = ;^(</'i + ^02)- The lattice 
action then reads 

ri A ^ 



.,x+i 



i=l 



cosh fl <pa,x4>a,x+i + ^ ^^^^ 1^ £ab<t>a,x4>b,^ 



X + 4: 



(2.3) 



We use the antisymmetric tensor Eat, with £12 
over repeated indices is implied throughout. 

Since the Boltzmann weight in the partition function, 



£21 = 1, £11 = £22 = 0, and summation 



Z 



D(j)iD(j)2 e 



(2.4) 



is complex, the theory has a sign problem and one cannot rely on importance sampling. 



Writing the weight as | e | e 



S\pi<4) 



e ^^e , one may consider the phase quenched theory 



^pq = J Dcl)^DcP2\e-^\ = J DcP^DcP2e-^^, 



(2.5) 



where is the real part of the action in Eq. ( p.3[ ), i.e. the term proportional to sinh/i 
is dropped. By analysing the average phase factor in the phase quenched theory, (e*'^)pq, 
it was shown in Ref. |l|] that this theory has a severe sign problem: at nonzero chemical 
potential the average phase factor goes to zero exponentially fast in the thermodynamic 
limit. 

We use stochastic quantization. The Langevin equations for the fields (pa read 

6S\ 



d 

-QQ(pa,x{9) 



S(t>aA^) 



+ Va,x{0), 



(2.6) 
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where 6 is the Langevin time. The noise r] is Gaussian and normahzed as 
Since the force in Eq. ( |2.6D is complex, the fields are complexified as 



ba^(p'^ + i(t>\ (a = 1,2). 



The complex Langevin equations we consider in this paper then read 



de 

The noise is chosen to be real. The drift terms are defined as 



K. 



K. 



R 
a,x 

I 



-Re 



5S 



Tm ■ 



SS 



(2.7) 
(2.8) 

(2.9) 
(2.10) 

(2.11) 
(2.12) 



and read explicitly 

= -[2d + m' + X {4>l^ - 41)] c^l, + 2\4>lA,x4>ix + E 



+ 



+ cosh fl 



+ 



a,x+4: ^a,x—4 



+ sinh^ Sab (4>\ 



x+'i ^b,x—i 



(2.13) 



ii 



+ 



a,x+i ^ a,x—i 



sinh/i Eab [4>f. 



-6^ ~ 

x+A ^b,x-4: 



(2.14) 



Observables are written in terms of the complexified fields ( |2.8D as well. We consider the 
square of the field modulus. 



1 



1 



|2 _ . i /-^Ra ^l2 

2^" 2 



+ i< 



(2.15) 



and the density (n) = {1/Q)dln Z/dfi, given by n = (1/^2) with 



n 



■X = {Sab Sinh (U - lEab cosh fl) (t>a,x4'b,x+i 



{6ab sinh ^ - lEab cosh /i) 



kR aR 



kl aJ 



a,x^b,x+4: ^a.,^^b,x+4: 



+ i 



kR jJ 



+ 



kl aR 



a,a;^fe,x+4 ' ^ci,^^b,x+4: 



(2.16) 



After complexification all observables have a real and imaginary part. 

Employing that the noise is random, we observe that the Langevin equations have the 
following symmetry. 



^l.xi 



'l,x 



^ix 



^ix, 



(2.17) 
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for all X and similar with 1 and 2 interchanged. Under this transformation, the drift terms 
change as 

^^,x ~^ ^^,x ~^ ^^^x-i ^l,x ~^ ^l,xi ^\,x ~^ ~^2,x- (2-18) 

Correlation functions odd under this transformation should vanish after noise averaging, 
which implies that 

{4>f,x(l>2',y) = {4'l,x4>l,y) = {(I>f,x4'l,y) = {4'2',x4'2,y) = 0- (2.19) 

The nonzero combinations are 



{€,x4>ly) 



^l,x4>l,y) 



t>a,x4>iy) 



(2.20) 



Applying this to the expectation values of the observables in Eqs. ( |2.15 , 2.16| ), we find that 
they are purely real. This is indeed what was observed numerically in Ref. M. 



3. Ignoring interactions 

In order to set the stage for the mean field analysis, we first solve the Langevin dynamics 
without interactions (A = 0), allowing for a detailed understanding of convergence and 
stability properties in the Silver Blaze regime. 

3.1 Standard results 

We start by summarizing the results obtained in the standard field theory approach (see 
e.g. Ref. pO]). After going to momentum space, according to 



'^a,x 

P 



^e'P^cPa,p, (3.1) 



where pi = 27rni/Ns, with —Ns/2 < rii < Ns/2, and p4 = lirn/^/Nr, with —Nt-/2 < n4 < 
Nt-/2, the action ( ^^ ) reads 



S = \4>a-p {5ah\ - EabBp) (f>b,p = ^ ]-(j)a-pMab,p(t>b,p, (3.2) 



p p 



where 



and^ 

3 

Ap = rr? + 4 sin^ ^ + 2(1 — cosh ^005^4) , -Bp = 2 sinh/x sinp4. (3.4) 

i=\ 



In the formal continuum limit Ap ^ m — fi + + p , Bp ^ 2fip4, 
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Note that A-p = Ap, B-p = —Bp, and that Mp is nonhermitian. The phase quenched 
theory is obtained by taking Bp = 0. Up to an irrelevant constant, the logarithm of the 
partition function is 

InZ = -ij^lndetMp = -^Y^H^ + ^p)- M 
p p 

The observables we are interested in are given by 

p p p 

and 

<"> = «lSr = -«5 + ■ 

where A' = dA/dfi = — 2 sinh/x cosp4, B' = dB/dfi = 2 cosh/xsinp4. 

As always, the severeness of the sign problem is estimated by the average phase factor 
in the phase quenched theory, given by the ratio of the partition functions of the full and 
phase quenched theories (^^, |2.5|), 



Z 



(e^^)pq = ^ = e-^^/, (3.8) 



where A/, the difference between the corresponding free energy densities, is given by 

A/ = -lln^ = ^yin^J±:S. (3.9) 



Note that this can be easily generalized to arbitrary powers of the phase factor in theories 
with nonvanishing phase factors.^ 

Finally, since the eigenvalues of Mp in the action (|3.2| ) are Ap±iBp, the theory without 
interactions exists provided that Ap > 0. This yields the standard stability criterium for a 
free Bose gas at finite chemical potential, 

Ap>0 ^ 4sinh2|<m2, (3.10) 

corresponding to //^ < m? in the formal continuum limit. We restrict the analysis below 
therefore to the case that Ap > 0; this is the Silver Blaze region. 



"Define the partition function Ze = J D(t>iD(t>2 le^^le""^. Then {e^"'^)i = Z„+i/Zi = exp{-i}Af„+i,e), 
with 

. , 1 , Z„+e 1 Al + (n + lfB^ 

A/„+,.^ = --ln^ = — 2^1n ^2+^2^2 ■ 
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3.2 Continuous Langevin dynamics 



We now solve the complex Langevin equations to compare the outcome with the results 
given above. The Langevin equations (^]^, 2.1C| ) read in momentum space 



d_ 
dd 



dd 



where the noise is normalized as 



Ignoring interactions, we find for the drift terms 

I 



K. 



-Ar. 



^a,p ^ ^^P^a,p •'^p^ab'Pb,pJ 

where Ap and Bp are defined in Eq. ( |3.4D above. In terms of 



(3.11) 
(3.12) 

(3.13) 



(3.14) 
(3.15) 



f Ap -iBp\ 
Ap iBp 



iBp Ap 







\-iBp Ap J 
the Langevin dynamics is written as 

do ^ 



-Mp<^p + Ep. 



( Vi,p\ 


m,p 
V J 



(3.16) 



(3.17) 



The matrix can be diagonalized by an orthogonal transformation and has doubly de- 
generate eigenvalues Xp = Apziz iBp. The solution of the Langevin equations is 



i>a,piO) = e-^"' [cos{Bpe)4>lp{0)+ism{Bpe)e,b4A^)] 

+ / dse-^"^'-'^ cos[Bp{e - s)]ria,p{s), 
Jo 

i>ipi^) = e"^^' [cos(i?p0)0i,p(O) -ismiBpe)eab4>lp{0)] 

-i / dse-^^'^^-'hinlBpie - s)]eabTlb,p{s), 
Jo 



(3.18) 



(3.19) 



R. I 

where (/>a,p(0) denote the initial conditions. 

We are now in a position to discuss the convergence properties of the Langevin process 
in the limit of large Langevin time. First we note that there is independence of initial 



conditions provided that Ap > 0, i.e. in the region of interest here. Taking (j)^^p{0) = 0, we 



find for the two-point functions, after using Eq. ( |3.13 ) and performing the Langevin time 
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2A,Al + Bl 



2Al + Bl 



pp 



2Ap + i?2 



integrals, 



«-pW<p'W> = '-j^^^i^Bp-e-'^^'[BpCos{2Bpe) + Apsin{2Bpe)]y (3.20) 

Most of the terms vanish in the hmit that 9 ^ oo, again provided that Ap > 0. The 
surviving terms are 



_g-2Ape ^^2 ^ ^2 _ ^2 (,os(25p0) + ApBp sm{2Bpe)] 



hm 



2^p ^2 -|- ^2 



hm (</.i,_p(0)</)[y(0)) ^ {^i-p4,p') = ^abS, ■ ^ 



^^2^p^2 + 52' 



hm ( 



6^ 



I I i B 

'b,p'i^)) = {4>a-p4>b,p') =£abSpp'- 



2AI + BI- 



(3.21) 



The structure of these two-point functions is in agreement with the symmetry ( 2.17 , 2.2C| ) 
discussed above. 

For the observables we find the following. The square of the field modulus ( ^.15 ) is 
given by 



-y 



An 



(3.22) 



which agrees with Eq. ( |3.6| ). After going to momentum space, the density (2.16) reads 
^"'^ ^ ((5„6sinh^cosp4 + Eab cosh ^ sin ^4) {(t>^_p(t>f^p - <fi-p(t>l^p + 2i4>^_p(t>l .p) 



E 



sinh ^ cos 



Aj + Bl 



cosh /i sin p4 



Br, 



(3.23) 



which agrees with Eq. (p.T]). Note that all two-point functions in Eq. ( 3.21j ) contribute to 
this answer. 

We conclude therefore that the Langevin process is independent of initial conditions 
and converges to the correct result in the limit of infinite Langevin time, provided that 
Ap > 0, as required in the Silver Blaze region. Moreover, the complexification is essential, 
as exemplified by the observables above. 
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3.3 Discretized Langevin dynamics 

We proceed by briefly considering the Langevin process after discretizing Langevin time as 
9 = ne, where e is the Langevin time step. The discretized Langevin equations are 

4>a,pin + 1) = + e<p(n) + V^r?„,p(n), (3.24) 

0i,p(n + 1) = cPi^n) + eKl^n), (3.25) 

and the noise obeys {r]a-p{n)r]i,y(n')) = 25nn'^ab^pp' ■ In the notation of Eq. ( p.l6| ), these 
equations are summarized as 

$p(n + !) = (!- eMp) ^p{n) + ^fe Hp(n), (3.26) 

and solved by 

n-l 

%{n) = (1 - eMpT ^p(O) + (1 - eMpf-''' Hp(i), (3.27) 

i=0 

where $p(0) is the initial condition. Convergence is determined by the eigenvalues of Mp. 
This yields the condition 

|l-eAp|<l, \p = Ap±iBp, (3.28) 

resulting in the constraint 

Ap-^-{Al + Bl)>Q. (3.29) 

We find that the convergence criterium is modified by an explicit stepsize dependence. 
However, this restriction is not special for the complex Langevin process |41, |3^. Consider 
real Langevin dynamics at zero /u or in the phase quenched theory. In both cases Bp = 
and the criterium reads 

0<^Ap< I. (3.30) 

Since Ap is maximal at the edge of the Brillouin zone {p = vr), this yields the modest 
constraint (for /i = 0) 

e < 73 n • (3.31) 

At nonzero chemical potential, this constraint is modified to 

2 

^ ^ 4(i + m2 + 2(cosh/x-l)' ^^"^^^ 

both in the full and the phase quenched theory. Li the Silver Blaze region, where /U is 
bounded, this leads to only a slightly stronger bound on e. However, since this bound is 
determined by ultraviolet modes at the scale of the lattice cutoff, it is likely that a similar 
constraint holds in the high-density phase as well. In the limit that 3> 1, exponentially 
small stepsizes would eventually be required. It should be noted, however, that for such 
large chemical potentials lattice artefacts are severe. 

The solution (l3^ ) can be used to study finite stepsize effects in two-point functions at 
infinite Langevin time. We come back to this below using a more elegant approach based 
on fixed points of the Langevin equations. 
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4. Fokker-Planck equation 



In order to better understand the Langevin process, we now study properties of the asso- 
ciated distributions. 

Consider first the Langevin process (|2.6D and the distribution P[(/>,0], defined via 



(4.1) 



where the brackets on the LHS denote noise averaging. The distribution satisfies the 
Fokker-Planck equation (in continuous Langevin time) 



{0[^,6]),= / Dct>P[<P,e]0[<P], 



dP[4), ^] _ ^ ^ 



As always, the index a = 1, 2 is summed over. The stationary solution 



(4.2) 



(4.3) 



always exists. However, since the action is complex, this is not the probability distribution 
for the complex Langevin process. 

More relevant for the complexified process (2.£, 2.101 ) we consider here, is the real 



distribution p[(jp' ,(f} ,6], defined via |26| 



(4.4) 



This distribution satisfies the extended Fokker-Planck equation 



de 



(4.5) 



If stochastic quantization is applicable for complex actions, the two expectation values 
dill ) and (U) should be equal 

We focus on the stationary solution of Eq. (|4.5|) and henceforth drop the 9 dependence. 
Ignoring again interactions, the stationary solution should satisfy 



E 



K. 



R 

a,p 



-K, 



0, 



(4.6) 



where the drift terms Ka,p were given in Eqs. (|1|, |3l|). Explicitly, this reads 

E 



s s s 



a,p 



R jJl 



0. 



(4.7) 
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Based on the structure of the equation, the solution can be written as 



R ^li 



A^exp 



E 



aR j,R 



a,p + Pp4>l,-p4>i,p + '^i£ab7p4>a,-p4>l,p) 



(4.8) 



where is a normahzation constant. Inserting this expression in Eq. (4.7) yields the 
coefficients 

A 



ap = Ap, f3p = -^{2Al + Bl) 

p 



Ip = -jr- 



(4.9) 



We have therefore found the stationary distribution corresponding to the complex Langevin 
process in the noninteracting case.^ Note that since 7_p = — 7p, the distribution is real in 
real space, as it should be. 

We now verify that this stationary solution is indeed the distribution corresponding 
to the Langevin process in the limit of infinite Langevin time. Performing the Gaussian 
integrals, we find the partition function 



1 



(4.10) 



where J\f is an irrelevant constant and 



A'. 



p 

The two-point functions that follow from this distribution are 



(4.11) 



^a,—pr'a,pl 



^a,—p4'a,p) 



p^b,p/ 



dlnZ 


Pp 




_ 1 2^2 + ^2 


dap 






Ap Ap + Bp 


dlnZ 


ap 




_ 1 


dPp 


C^pPp 




Ap A^ + B^ 


dlnZ 


-27p 


_ —'^Bp 


dip 


ap(3p 




A^ + B^ 



(4.12) 



These agree exactly with the results obtained by solving the Langevin equation, cf. Eq. 

The theory with the probability distribution p[(j)^,(l)^] exists provided that the eigen- 
values of the quadratic form in Eq. (|4.8|) are positive. We find the eigenvalues to be 



A. 



i (ap + pp± ^ {ap - f3pf + A^l) 
Al + Bj[^Al + Bl±Ap 



Ap 

p 



(4.13) 



These are positive provided that Ap > 0. The criterium that determines the convergence 
of the Langevin dynamics also emerges in the stationary solution of the extended Fokker- 
Planck equation, as expected. 



^See Refs. ^ for other examples. 
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Let us discuss some more properties of the distribution ( [4. 81) . First we note that the 
distribution is highly nonlocal in real space and does not allow for e.g. a derivative expan- 
sion, due to the division by Bp = 2sinh^sinp4 in the coefficients ([4. 9]). We find therefore 
that the complexity of the original local weight has been traded for the nonlocality of 
the real probability distribution. However, this nonlocal behaviour is expected: it follows 
from the Langevin equations that the modes with P4 = are purely real, i.e. (/)^ (p4=o p) ~ 
This is enforced in the probability distribution by the singular behaviour as p4 — > 0. For the 
same reason the limit — > is singular, since there is no need to complexify the dynamics 
in this case and the distribution for the (j)^ modes should reduce to a delta function, 5(</>^). 
These considerations fix the dependence on Bp. 

In conclusion, we have found the stationary solution of the extended Fokker-Planck 
distribution. The real distribution is nonlocal and singular in the limit that fi,P4 0. 



5. Mean field approximation 

We now return to the interacting theory, with discretized Langevin time 6 = ne, and 
consider the two-point functions 

GX(^) = «_p(n)4p(n)), 

Gilpin) = {4'l-p{n)cP{p{n)). (5.1) 

Using the Langevin equations ( |3.24| , |3.25| ), we find that these correlation functions evolve 
according to 

GXpi^ + 1) = GXpin) + e{<t>l.p{n)Kf;p{n) + <_p(n)0^p(n)) 

+^\Kf^-p{n)Kf;p{n)) + e(r?a,-p%,p), 
Gab,piri + 1) = G^,^p{n) + e{4>l_p{n)Klp{n) + Kl.p{n)4>{p{n)) 

+e\Kl_p{n)Kl^{n)), 
G^lp{n + 1) = Gilpin) + e«_p(n)i^,\p(n) + Kl_p{n)4>lp{n)) 

+e\Kl_p{n)Klp{n)). (5.2) 

Here we used that {r]a-p{n)<jx^'^{n)) = 0, since the fields at time n do not depend on the 
noise at time n. The terms proportional to are finite stepsize corrections. We then look 
for fixed points of the Langevin equations, and put 

GX{n + l) = G'^^p{n), (5.3) 

etc. This yields the fixed point equations 

{(t^l-pKlp + i^^p</>^p) + e{Kl_^Kl;) = -26at, 

{4>i-pKlp + <-pAp) + '(K-pKp) = 0' 

+ Kf^-p4,p) + <^{Kl.pKlp) = 0. (5.4) 
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To implement a mean field approximation and find a self-consistent set for the two- 
point functions ( |5.1| ), we factorize the interaction terms. Consider for example the term 
^f'x^b x^^x appearing in the drift term ( 2.14| ). We write 

4>f,x4,A ^ + i^b,x4>lxHl,x + {4>l,x€,x)^t- (5.5) 

Using the notation 



G^^in) = «,.(n)<,(n)) = 1 J] GS,'^(n), 

p 

Glin) = {4>ix{n)4>lM) = ^E^aip(^)' 

p 

G^lin) = {^lM4,xin)) = i ^ G^in), (5.6) 



p 



the drift terms in the mean field approximation read 

Kf,P = -[Ar + X {GT - Gil)] ^l, - 2A {G^ - G^) 

+ [iBpS,, + 2A (GSJ + Gf^)] 4>{p + 2\GWa,p, (5-7) 

- [iSpe^fe + 2A (gSI + G^i)] C - 2AG^,V^p. (5.8) 

These can be further simplified by noting that both G^^{n) and G^^fj{n) are proportional 
to 6ab, for all Langevin times. We write therefore 

G^^in) = 6,,G^^{n), 

Glin) = 5abG'\n), 

GSI(n) = 5abG^\n) + eabG^\n), (5.9) 
such that the drift terms reduce to 

K^,P = -Ma,p + iBpeab4>{p + C4>ip, (5.10) 

Kl,p = -■^p4>l,p - iBpeab(t)f;p - C(l)a,p, (5.11) 

with 

Ap = Ap + 4X {G^'^ - G^^) , C = 8XG^\ (5.12) 

Since Ap{n) and C{n) depend explicitly on the Langevin time, the time-dependent mean 
field Langevin equations cannot be solved analytically.^ 

We therefore look for fixed points. After substituting Eqs. ( 5.10| , 5.11) in the fixed 
point equations (|5.4| ) and performing some algebra, we find that at the fixed point C = 
and that the two-point functions can be decomposed as 

^ab,p — Oab^p ; '-^afe,p — "ab^p j ^ab,p — ^^ab'^p ■ {0.16} 



®In fact, the problem is now very similar to that of nonequilibrium field dynamics using a self-consistent 
mean field approximation in the equal-time formalism [M. 
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This is in agreement with the symmetry ( 2.17 , 2.20| ). The three fixed point equations (| 
then become 



•^pGp + BpGp ^ (^J[.pGp -\- BpGp -\- IjApBpGp ) — 1, 
•^pGp BpGp — (^jApGp + BpGp 2jApBpGp ) — 0, 



ApGf - \Bp [G^^ - G^i] -U[Al- Bl] Gf - ApBp [G«« - G^] ) = 0. (5.14) 



The solution is 



1 



G: 



Ap 2^{Ap + Bp^ 
RR qH — Ap{l- ^eAp) + ^eB. 



Al + Bl (1_ 1,^^)2 +1,2^2 ' 



G. 



RI 



1 



Bn 



1 - eAr, 



2Al + Bj^l_i,^^f + l,2B2- 
For vanishing Langevin stepsize this simphfies to 



rfRR I nil 

Up + (jTp 



1 



G 



RR 



G: 



G: 



RI 



1 



2^2 + ^2' 



(5.15) 



(5.16) 



ip - -p ' 

while in the phase quenched theory [Bp = 0) where real Langevin dynamics is applicable, 
the solution reduces to 

1 1 



G: 



RR 



Ap 1 



UA ' 

2t./Hp 



G" = G; 



RI 



0. 



(5.17) 



We find finite stepsize corrections linear in e, as expected [41]. Furthermore, we note that 
for large stepsize the denominator in the first line of Eq. ( 5.15 ) can go negative. However, 
this occurs precisely when the stability criterium ( 3.29 ) is violated (after the replacement 
Ap — > Ap) and is therefore excluded. 

The expressions in Eq. ( ^.16| ) agree precisely with the solutions ( |3.21| ) obtained by 
solving the Langevin equations without interactions, after making the mean field replace- 
ment Ap — > Ap. This replacement corresponds to the standard mean field approximation 
in which the mass parameter receives a tadpole correction, 



(5.18) 



or, in the notation of this section. 



m 



TT? + 4A {G^^ - G") 



4A 



Ar, 



Ad, + B^ 



(5.19) 



with 



Ap = M'^ sin^ ^ + 2(1- cosh fi cosp4) 



(5.20) 
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Figure 1: The lines represent the mean field results for in the full (left) and phase quenched 

(right) theories for various lattice sizes, taking m = A = 1. The vertical dotted line indicates the 
mean field estimate for the critical chemical potential. The data points are obtained with Langevin 
simulations jl]. 




Figure 2: As in Fig. |i| for the density (n). 



These equations define a self-consistent gap equation for M^. Given m and A, the gap 
equation can be solved numerically after specifying the lattice size. For example, taking 
m = A = 1, we find = 1.47 and = 0.119 at // = on a lattice of size lO'^. The 

critical chemical potential then follows from Ao = (or = 4sinh^ ^) and is found to 
be He = 1.15.^ 

We have solved the gap equation in the Silver Blaze region and used the outcome to 
compute {\4'\'^) and (n) in the mean field approximation as a function of chemical potential 
for different lattice sizes. The results are shown in Figs. |^ and |2| respectively, for m = 
A = 1. The vertical dotted lines indicate the mean field estimate of the critical chemical 
potential. In the full theory (figures on the left) the expected /i independence emerges in 
the thermodynamic limit. In the phase quenched theory (figures on the right), observables 
depend on fi, since there is no Silver Blaze feature, see Appendix Also shown in these 



Wlien X^O, fic follows from = (or = 4sinh^ |), yielding = 0.962 for m = 1. 



Figure 3: Left: as in Fig. for the average pliase factor in the phase quenched theory (c*'^)pq- 
Right: difference A/ between the free energy densities of the full and the phase quenched theories 
in the mean field approximation. 

plots are data points obtained from the numerical solution of the Langevin process, with 
stepsize e = 5 x 10~^ We observe surprisingly good agreement between the mean 
field and the nonperturbative results for all values of the chemical potential and all lattice 
sizes considered, indicating that the mean field approximation captures the most relevant 
interactions. 

In Fig. ^ we show the average phase factor in the phase quenched theory (e*'^)pq 
(left) and the difference A/ between the free energy densities, given in Eq. (|3.9| ), again 
after the replacement Ap — > Ap. As already mentioned, the sign problem is severe in the 
thermodynamic limit: taking e.g. A/ ~ 0.02 and a lattice volume Q = W^, we find that 
the average phase factor {e^'^)pq = e~^^^ is indeed exponentially small. 

To conclude this section, we note that it is straightforward to adapt the stationary 
solution of the extended Fokker-Planck equation, constructed in Sec. ^, to the mean field 
approximation discussed here. Since in the mean field approximation only two-point func- 
tions appear, the mean field probability distribution remains of the form (|4.8| ) with the 
simple replacement Ap — > Ap (or tv? — > M^). Existence of the Fokker-Planck distribution 
in the Silver Blaze region now requires Ap > 0. 

6. Conclusion 

In order to further understand the applicability of complex Langevin dynamics for theories 
with a complex action due to finite chemical potential, we have studied the relativistic 
Bose gas in the Silver Blaze region analytically. Ignoring interactions, we have investigated 
convergence and stability, and constructed the stationary solution of the extended Fokker- 
Planck equation. We explained why this real probability distribution is nonlocal in real 
space. Subsequently, interactions were included on the mean field level and the fixed 
point of the mean field Langevin equations with finite stepsize was derived. We gave a 
comparison between the mean field predictions and the nonperturbative numerical data 
from Ref. in the Silver Blaze region. Surprisingly good agreement was found for all 



values of the chemical potential considered, including finite size effects, indicating that the 
mean field approximation captures the most important effects of the interactions. We have 
demonstrated analytically that the sign problem is severe for lattice volumes used in this 
study. From the combination of results obtained here and in Ref. |T|, it can be argued that 
complex Langevin dynamics in the Silver Blaze region is well understood in this theory. 

One obvious next step is to extend the analysis to the high-density phase, which 
requires the introduction of the mean field ((^^^.) (note that (i;^^^^,) = 0). Finally, it would 
be interesting to apply mean field approximations to other theories as well, especially in 
combination with numerical studies. In particular, this would be useful for QCD with 



static quarks pjq| . 
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A. Dispersion relation 

The propagator corresponding to the action ( |3.2D is 

^ _ ^ah-^p + ggb-Bp /' A 1 ^ 

Dispersion relations follow from the poles of the propagator, taking p4 = iE-p. We find 



1 o\ /I 

cosh Ep[fj,) =11 + 

where 



cosh Ep{fj,) = ( 1 + -Ljp j cosh fi± y 1 + -tJp sinh /i, (A.2) 



'2 2 I A • 2 P' 



p 



4^sin2|. (A.3) 



This can be written as 



cosh Ep{^i) = cosh [Ep{0) ± fi] , (A.4) 

such that the (positive energy) solutions are 

Epifi) = Ep{0) ± fj,, (A.5) 

just as in the continuum theory. Lattice discretization effects only appear in the disper- 
sion relation at zero chemical potential, £'p(0). The critical fj, value is Hc = Eo{0) = 
2asinh(?TT,/2), so that one mode becomes exactly massless at the transition. 

The phase quenched theory corresponds to Bp = 0. In the formal continuum limit, the 
phase quenched theory is a theory with a real action and mass parameter — /i^. The 
dispersion relation is 

"*^-<''> = ^i;('+rO' 

corresponding to Epi^n) = m? — ^u^ + in the continuum limit, as anticipated. 

These results are easily extended to the self-consistent mean field approximation, where 
the mass parameter receives a tadpole correction and is replaced by = m? + 4A 
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